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Abstract. We investigate the scattering cross section of aligned dipolar molecules 
in low-temperature gases. Over a wide range of collision energies relevant to 
contemporary experiments, the cross section declines in inverse proportion to the 
collision speed, and is given nearly exactly by a simple semiclassical formula. At 
yet lower energies, the cross section becomes independent of energy, and is reproduced 
within the Born approximation to within corrections due to the s-wave scattering 
length. While these behaviors are universal for all polar molecules, nevertheless 
interesting deviations from universality are expected to occur in the intermediate 
energy range. 
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1. Introduction 

The energy dependence of scattering cross sections for atoms at ultralow collision 
energies is very simple. Either the cross section is nearly independent of energy, 
for distinguishable particles or identical bosons; or else the cross section vanishes 
altogether, for identical fermions. This behavior emerges in the limit where the 
deBroglie wavelength exceeds any natural length scale of the interatomic interaction, 
and the scattering is characterized by a single quantity, the s-wave scattering length a 
(alternatively, the p-wave scattering volume V p for fermions). Although a dominates 
the threshold scattering, nevertheless its value is not immediately obvious from the 
interaction potential, and must be determined painstakingly from experiments. 

By contrast, for low-energy collisions between polarized dipolar molecules, the near- 
threshold scattering is often approximately determined directly by parameters of the 
interaction potential. The interaction between two molecules of reduced mass M and 
dipole moment /i is characterized by a dipole length, given by D = Mfi 2 /h 2 . H Under a 
wide variety of circumstances, to be discussed in this paper, D is the dominant length 
scale and sets the threshold cross section, i.e., a ~ D 2 . Moreover, this circumstance 
holds for identical fermions as well as bosons. In this sense, the scattering of two dipoles 
is nearly universal at threshold, apart from possible modifications arising from s-wave 
scattering. 

For realistic collisions in present-day experiments, however, collision energies are 
not always in this threshold region. When the threshold region is left behind, there is 
significant numerical evidence to suggest that a universal behavior still emerges, and that 
the cross section scales as a ~ D/K, where K is the wave number of the relative motion 
PQ. The switch between the two types of behavior corresponds roughly to the natural 
energy scale of the dipolar interaction, Ed = fx 2 /D 3 = fr 6 /M 3 /z 4 . Below this energy 
nonzero partial waves contribute only perturbatively, and only at large intermolecular 
separation; whereas at higher energies, many partial waves contribute and the scattering 
is semiclassical. 

These two behaviors may thus be said to define the operational limit between 
"cold" and "ultracold" regimes of scattering for dipoles. We propose that the onset of 
cold collisions occurs when the temperature of the gas passes well below the molecule's 
rotational constant B e (or else its A-doublet splitting) so that the orientational degrees 
of freedom freeze out. These temperatures are typically in the mK-K range. This cold 
collision regime, in which semiclassical scattering occurs, persists until the temperature 
gets as low as Ed (typically nK-/iK temperatures for dipolar molecules). Temperatures 
below Ed define the ultracold regime, where true threshold scattering is apparent. 
Actual values of these temperatures depend strongly on the species considered and 
on the applied electric field. 

\ D is determined, roughly, by equating a typical centrifugal energy, h 2 /MD 2 , to a typical dipolar 
energy, n 2 /D 3 . This is the same reasoning that leads to the definition of the Bohr radius, by equating 
centrifugal and Coulomb energies for hydrogen. 
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Figure 1. The total scattering cross section a, averaged over all incident directions, 
for two distinguishable polarized dipoles. In the low-energy limit, a reduces to the 
Born approximation result, while at high energies it is given by a semiclassical eikonal 
approximation. Note, however, that boundary conditions were carefully chosen in this 
calculations so that the s-wave scattering length vanishes. 

Our goal in this article is to make these ideas precise. We will illustrate the universal 
behavior of dipole-dipole scattering at low temperatures in the two regimes, and, more 
importantly, we will see the circumstances under which this universality fails. Figure 1 
shows the basic elements of universality for dipole-dipole scattering. This figure plots 
the total scattering cross section a, averaged over all incident directions, versus collision 
energy. Both quantities are presented in terms of the "natural" units given above. The 
black curve is a complete numerical close-coupling calculation. At low energies, E < En 1 
a approaches a constant value that is well-approximated using the Born approximation 
(blue). At higher energies E > En, & falls off as 1/yE, and is given by the semiclassical 
eikonal approximation (red). Both approximations will be derived below. 

In this article, we present alternative close-coupling calculations that reveal 
deviations from this universal behavior. For example, at threshold, the cross section 
can deviate substantially from the Born result. This occurs when threshold resonances 
in s-wave scattering dramatically increase the cross section above the D 2 value (2J [3J 
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HI More interestingly, in the intermediate energy regime E m where neither 
approximation holds, our close-coupling calculations show that the behavior of cross 
section versus energy depends on details of the physics on length scales small compared 
to D. For this reason, it is conceivable that elastic scattering experiments that probe 
deviations from universality will be an important tool in unraveling information on the 
short-range physics of close encounters. Finally, we also discuss the angle-dependence 
of cold collisions, in which the incident direction of the collision partners is varied with 
respect to the polarization axis; this is a measure of the importance of the anisotropy 
of the dipole-dipole interaction. These results should serve as simple and accurate 
guidelines to low-energy collision cross sections needed to understand scattering or 
thermal equilibration in contemporary experiments. 



2. Formulation of the problem 

For a pair of dipoles with reduced mass M and polarized in the z direction by an external 
field, the two-body Schrodinger equation reads 

h 2 „ , , , ,1 — 3 cos 9 



V 2 — + V SR 



2M v R 3 



V> = Eip. (1) 



Here R = (R,9,<p) is the relative displacement of the two dipoles, and (//i) and (/x 2 ) 
are their induced dipole moments. These dipoles are field dependent, and their values 
are set by the field and the internal structure of the dipole. In the present context, we 
will take them as fixed quantities, and observe their influence on the scattering. By 
ignoring internal molecular structure, we are in effect modeling molecules trapped in 
their absolute ground state, as in an optical dipole trap. In particular, we assume that 
inelastic scattering is disallowed, a topic taken up in the next article of this issue [Sj. 

The potential term Vsr represents all the short-range physics, i.e., the potential 
energy surface of chemical significance when the dipoles are close together. This detail 
will be specific to each pair of collision partners considered. To simplify the discussion, 
we will replace the complex details of Vsr by imposing a boundary condition at a fixed 
interparticle separation Rq. We will assert that the wave function if> vanishes uniformly 
on this boundary, although alternative boundary conditions are certainly possible and 
indeed desirable [21 [3] . Having made this approximation, we ignore Vsr in the following. 

The resulting Shrodinger equation then admits a natural length scale D = 
MOuxXa^)/^ 2 , an d a natural energy scale E D = h 6 / M 3 (fi^ 2 (/j, 2 ) 2 ■ By recasting (Tj[|) 
in the scaled units r = Rf D, e = E/Ed, and by ignoring Vsr, we arrive at the universal 
Schrodinger equation (assuming that the molecules have the same orientation relative 
to the field axis) 

1 2 2C*20 



V 



, v j V = e0, (2) 

where C2o{9,4>) = (3cos 2 # — l)/2 is the standard reduced spherical harmonic [7j. 
Although the equation has a universal form, nevertheless its solutions may depend 
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on details of the short-range physics, represented here by the cutoff radius r = Rq/D. 
Our main objective is to explore circumstances under which solutions are universal, i.e., 
independent of tq, under the conditions of modern cold molecule experiments. 

As a preliminary argument in this direction, let us consider the relative sizes of 
the two characteristic lengths D and Ft®. For most atomic and molecular species, the 
size scale Rq below which short-range physics can matter is on the order of the van 
der Waals length, which is typically ~ 100 a , which is also the distance at which the 
internal fields generated by the dipoles are comparable to the applied field (at least 
for typical laboratory field strengths). By contrast, for molecules with typical 1 Debye 
dipole moments, D will be orders of magnitude larger. For example, in the OH radical, 
with fi = 1.68 Debye, -Dqh = 6770 a ; for the representative alkali dimer KRb, with 
/i = 0.566 Debye, D-^Rb = 5740 ao; and for the highly polar LiCs molecule with /i = 5.5 
Debye, we find -Dlics = 6 x 10 5 ao- Correspondingly, the characteristic energies are 
low: E D £,n = 445 nK, while E DKRh = 83 nK, and -E^LiCs = 7 pK. In making these 
estimates, we assume the molecules are completely polarized; D gets shorter, and En 
higher, if they are only partially polarized. |§| 

The finite size of Rq can destroy universality in three ways. First, at low energies, 
we will see that scattering cross sections are of order D 2 . This would be spoiled by 
the geometric cross section oc -Rq, if Rq were comparable to, or larger than, the dipole 
length D itself. However, we have just argued that this is usually not the case for polar 
molecules. Second, there is the possibility that the s-wave scattering length a is larger 
than D, and this would also alter the universal result at low energy. Indeed, this very 
idea has been invoked as a means of tuning the interaction between dipolar molecules 
[TT1 [T2l [TBI [T4l [15] or else as a tool for probing details of these interactions [2j. Often, 
the effect of the scattering length is non-negligible, as we will see below. 

The third instance where Ro may matter is in the extreme high energy limit, where 
the universal cross section falls as D/K, where K is the wave number. In this case the 
cross sections will again tend to the geometrical ~ -Rq for our artificially imposed hard 
wall, and will dominate at energies where -Rq > D/K, which translates to about 4.5 
Kelvin in OH, well above Stark decelerator energies, and at which point other degrees of 
freedom of the molecule are relevant. Thus non-universal behavior may not be a concern 
at "high" energies either, from the standpoint of current experimental investigations. 

In this article we compute scattering cross sections for dipoles over a wide range of 
collision energies. We do this in three ways: 1) a fully numerical close-coupling expansion 
of the wave function in partial waves; 2) a Born approximation that exhibits the correct 
universal behavior (a ~ D 2 ) in the ultracold limit (E < Ed)', and 3) a semiclassical 
eikonal approximation that exhibits the correct universal behavior (a ~ D/K) in the 
cold collision regime (E > E^,). We briefly describe these methods in the following 

§ For species whose magnetic dipole moment mediates the interaction, these scales can be quite 
different. For atomic chromium, for instance, one finds Dcr = 27 ao, far smaller than its natural 
scattering length. A substantial body of literature now treats the complete details of the Cr-Cr 
interaction potential, for example [HI 03 QH] • 
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subsections. 



2.1. Close- coupling formalism 

We expand the total wave function into partial waves in the conventional way: 

^(r,e,<f>) = -J2Fi m (r)Y lm (9,(f>), (3) 

Ira 

where the radial functions F\ m satisfy a set of coupled-channel Schrodinger equations 
1 d 2 F lm 1(1 + 1) 2 (m) 

~2l^ 1^ ,m 7y "' I'Vrn-el'lm, (A) 

and the coupling matrix element is given by [7] 

= (lm\C 20 \l'm) (5) 

Z 2 




= (-1)"V(2Z + 1)(2Z' + 1) , 

Owing to the cylindrical symmetry of the Hamiltonian, the angular momentum 
projection m is a good quantum number, and we can solve a separate set of coupled 
equations for each value of m. However, the boundary conditions of the wave function 
may not respect this symmetry, i.e., the incident wave could arrive from any direction, 
not just along the symmetry (z) axis. Thus we consider a complete sum over m in the 
wave function (jHJ). Similarly, symmetries of the 3-j symbols in (jSJ) guarantee that each 
partial wave Z is coupled only to the partial waves V = I, I ± 2 by the dipole interaction. 
Thus we can consider even partial waves separately from odd partial waves, and will do 
so in the following. For identical particles in the same internal state, these situations 
correspond to bosons and fermions, respectively. 

The equations (J4j) admit as many linearly independent solutions as there are 
channels (Zm). Individually, they are defined by the boundary conditions (for each 
m) 

Ft(r = r ) =0 

F}^(r - oo) = 5 w e-^ kr - 1 '^ - S$?> e^- 1 '*™ , (6) 



where k = y2e = DK is the wave number in dipole units. These scattering boundary 
conditions serve to define the scattering matrix . From this matrix one can construct 
the scattering amplitude describing scattering in direction kf — (6,<j)) from an incident 
direction ki [T6| [T7]: 

f(k k f ) = ~E i%m(h)T^ ] r l % m (k f ), (7) 

K IV m 

in terms of the T matrix, T = i(S — I). By integrating over the final directions, we 
arrive at the total cross section for dipoles incident along kf 

0tot(&i) 



D 2 



d(j)d(cos6)\f\ 2 

An - - 

T mki,ki). (8) 
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This last line is the familiar optical theorem result. This is the type of cross section that 
can be measured in cold beam experiments, where, say, one species is trapped and the 
other is incident on the trap from the terminus of a Stark decelerator. If the molecules 
are magnetically trapped, then an electric field can be applied an at arbitrary angle 
relative to the collision axis |18j . 

Finally, the total cross section, integrated over an assumed isotropic distribution of 
initial directions, is 

D 2 J D 2 

^ ll'm 

This cross section is more relevant to in situ collisions in a trap, which serve to re- 
thermalize the gas and provide evaporative cooling. 

Numerical solutions to the coupled-channel equations (j3J) are determined using a 
variable stepsize version of Johnson's algorithm [20J. To ensure convergence of total 
cross sections, we include partial waves up to I ~ 100 at the highest collision energies of 
e = 10 4 . Vice versa, at the lowest collision energies we can get away with partial waves 
up to / ~ 30, but must apply the boundary conditions (J6j) as far out as r = 20, 000. 



2.2. Born approximation 

At the lowest collision energies, the Wigner threshold laws are well known to be different 
for dipolar interactions than for, say, van der Waals interactions. The elastic scattering 
phase shift 8\ in partial wave I > 0, due to a potential with l/r s long-range behavior, is 

m 

tan 5, ~ Ak 2l+1 + Bk s ~ 2 , (10) 

for some constants A and B that depend on short-range details. The first term in ffTUl) 
arises from the action of the short range potential, while the second is due to purely 
long-range scattering outside the centrifugal barrier. For the van der Waals potential 
(s = 6), both contributions go to zero faster than ~ k at zero energy, and do not 
contribute to the threshold cross section. However, for the dipole-dipole interaction 
(s = 3), the second term is ~ k for all partial waves. The contribution to the cross 
section, cx siia 2 5i/k 2 , is then independent of energy in all partial waves, and this cross 
section arises from long-range scattering. 

This circumstance leads to the applicability of the Born approximation in threshold 
scattering of dipoles [21], [13], [H] [221 [23] EES]. At ever lower collision energy, scattering 
occurs at ever larger values of r outside the barrier. But at long range the dipole-dipole 
interaction cx 1/r 3 is weak, and the perturbative Born approximation is applicable. 
These remarks do not apply, however, to s-wave scattering, where there is no barrier. 

The Born approximation for the scattering amplitude reads 

f{k h k f ) = ~l d 3 re- rk r r y d (r)e rk ^. (11) 

Z7T J 
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Replacing each plane wave by its standard partial wave expansion and re-arranging 
yields 

f(h, k f ) = -^1 d ' r ("4) c *><f) 47r E r e Y Vm 4 f )Y? m ,(?)Mkr) 

I'm' 

x^J2i% m (k)Y lm (r)ji(kr). (12) 

lm 

Consolidating the integrals into radial and angular varieties, we arrive at 

f(k u k f ) = 2tt £ i%* m {k)C\fY n a- v Y Vml {k f ), (13) 

ll'm 

where is the angular integral fl5J) defined above, and r M / is the radial integral 

v ' 2,Ji(kr)ji'(kr) 



Tw = 8 r'alr 



tiT((Z + Z')/2) 



r((-z + /' + 3)/2)r((/ + z' + 4)/2)r((z - Z' + 3)/2) ' 
Z' = z 

32 Z' = Z + 2 



3(Z+l)(/+2) ' 

Comparing the Born result f|T3|) with the expression (jTJ) (the factor of 8 is intended to 
simplify this) identifies the T-matrix in the Born approximation as 

T M,Born = _ kC M Tll ,_ (15) 

The Born approximation must be applied with a caveat. For purely s-wave 
scattering, where Z = I' = 0, the matrix element Cqq vanishes, and so therefore 
does the T-matrix element T q ' Born . The Born approximation is therefore mute on 
the question of s-wave scattering. As argued above, s-wave scattering, described by 
a scattering length a, is part of the non-universal behavior of scattering anyway. To 
produce realistic scattering results, it is possible to supplement T Born with an empirical 
s-wave contribution, which is determined from the full close-coupling calculations 

[IB EH USE]. 

Summarizing all these results, the threshold cross section in the Born 
approximation, averaged over incident directions, is given by the incoherent sum 

^Born 71 i/T-,(m),Borni2 f-i f?\ 

Win 

Evaluating the sums, the even and odd partial wave cross sections at threshold can be 
given as 

(j^ orn = 1.117 D 2 + 4vra 2 

o-° Born = 3.351D 2 , (17) 

where Ana 2 allows for the existence of a scattering length a, which is not determined 
in the Born approximation. For identical particles, the cross sections ( 1T7|) must be 
multiplied by 2, as usual; for distinguishable particles, both even and odd partial waves 
are possible, and these cross sections are to be added. 
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The low-energy limits (1171) were independently verified in near-threshold 
calculations using a coupled-channel adiabatic representation [U [5], which also 
illustrated how threshold angular distributions are affected by the competetion between 
long-range (Born) coupling and s-wave scattering resonances, as discussed in Section 4. 



2.3. Eikonal approximation 

At sufficiently high energies, a semi-classical analysis yields the simple D/K scaling of 
the cross section. Note, in reference to Fig. (1), that the deBroglie wavelength becomes 
smaller than the natural dipole length scale, 2ir/K < D, when E > 2tt 2 Ee>. This 
marks the onset of semi-classical scattering. This semi-classical onset can lie at fiK 
temperatures, or even colder, owing to the large dipole length scale and therefore the 
small value of E^. Many partial waves contribute to the scattering amplitude in the 
semi-classical regime, and differential cross-sections are increasingly concentrated in the 
forward direction. 

The eikonal method was long ago developed to find approximate scattering solutions 
of wave equations such as ([2]) valid in the semi-classical or ray-optics limit, in which the 
potential is assumed to vary little on the scale of the wavelength. A derivation of the 
eikonal wavefunction will not be given here, as it can be found in familiar texts [21] 
and in a comprehensive review article by Glauber [25]. Suffice it to say that a phase- 
amplitude ansatz, coupled with the assumption of a slowly-varying amplitude, leads 
directly to an approximate wavefunction 



ip(r) — e i ' r exp 



V(b, <j>, z')dz' 



(18) 



i 
k 

Following Glauber, this wavefunction is expressed in cylindrical coordinates with a new 
quantization (or z) axis aligned with the average collision momentum, fc avg = (ki + kf)/2. 
The cylindrical radius about this axis, b, can be associated with a classical impact 
parameter: is the azimuthal angle about the quantization axis. Due to the shift of 
quantization axis away from the direction of the applied field, S, we now note the field 
direction explicitly in the potential V(r) = [1 — 3(f ■ £) 2 }/r 3 . 

The two factors in the eikonal wavefunction are familiar in the context of the 
one-dimensional WKB method, where they coincide with an expansion of the WKB 



phase i J dz y2(e — V) ~ ikz — (i/k) J dz V to first order in V/e. Accordingly, 
we anticipate that the eikonal method will be most accurate when the incident energy 
is large compared to the magnitude of the dipole-dipole interaction, a more stringent 
criterion than the semi-classical constraint noted above. 

The analysis of scattering amplitudes associated with the approximate eikonal 
wavefunction is simplified by a judicious choice of coordinates. Note, in particular, that 
the momentum transfer q = ki — kf is orthogonal to the quantization axis defined by 
/c avg . For simplicity, we define an x-axis along q, in which case the y-axis is orthogonal to 
the collision plane and lies along k avg x q. In this reference frame, the impact parameter 
is written in vector form as b = 6cos(0)x + 6sin(0)y, and the relative displacement of 
the dipoles is r = b + zk avg . 
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Insertion of the eikonal wavefunction, valid where the potential is non- negligible, 
into the integral equation for scattering leads to the eikonal scattering amplitude 



f Ei (kt. fa = — 
J K f ' 11 2ni 



bdb 



iqbcos(<f>) 



- 1 



(19) 



where k = \kj 



X(b) 



k 



2e, q — \q\ — 2fcsin(# s /2), and where the eikonal phase is 
dz' 1/(6,0,2') 



With the explicit form of the dipole-dipole potential 



V(b,<p,z) 



, (b ■ S + zk avg ■ £f 
b 2 + z 2 



(&2 + ^2)3/2 

the phase is readily evaluated; setting a = z/b one finds 



(20) 



(21) 



1 



X= - 



da 



kb 2 
- 3(b-S) 2 

~~ 3(/c av g • S) 



-oo {l + a 2 f/ 2 
da 



1 + a 2 ) 5 / 2 
50 a 2 da 



'1 + a 2 ) 5 / 2 



The integrals are straightforward, giving 

2 



X 



Referring the electric field to our coordinate axes (x — q,y — k. 

£ = sin a cos (3x + sin a sin (3y + cos ak avg 
the phase is simply 



avg 



x q, k, 



avg; 



X(b,<f>) 



kb 2 



sin 2 a cos(20 — 2(3) 



The eikonal amplitude now has the form 



2m 



b db 



iqb cos <j> 



exp { i—— sin 2 a cos(20 — 2(3) 
kb 2 



(22) 

(23) 
(24) 
(25) 
(26) 



Explicit evaluation of the resulting integrals has proven quite difficult. However, 
to extract total cross-sections, Glauber's general proof of unitarity of the eikonal 
approximation [25] permits use of the optical theorem, (jBJ). 

For forward scattering, q = and the first phase vanishes identically. Expressing 
the result in terms of the orbital angular momentum I = kb gives 

1 



/ ki) 



2i\ik 



I 2k 

exp < %— sin 2 a cos(20 — 2(3) 



(27) 



but note that this appears undetermined since q and, consequently, x and y are 
not defined when q = 0! In this limit, a = arccos(A; • £) is well-defined, but 
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(3 = arctan(y • S/x ■ £) is not. Fortunately, it is easy to show that the azimuthal 
integral is independent of /3, with the result 



/ % (k h h) - — 



Jo -zr sin 



a 



From the optical theorem, the total cross section is then 



a 

15 2 



poo 




/ I dl 


1 - J o (jZ sin2 a ) 


Jo 





(28) 



(29) 



This result provides some insight into partial wave analysis in the semi-classical regime, 
which approximately separates into two regions: For I < v / ^sin(a), the integrand is 
nearly linear in /, while for larger / it declines steeply as A; 2 sin (a) /I 3 . Using these 
approximations, the integral evaluates to k sin 2 (a), the exact result given below. 



More carefully, we set 
2k 



sin 2 a 



d£ 



-k sin a 



ds 



P 8' 

and so express the total eikonal cross section in terms of a dimensionless integral 
rLi 47rsin 2 a r°° ds ri T , 



a 







The integral evaluates to unity, with the result 

k 



a. 



totV 

D 2 



ki ■ S 



(30) 



(31) 



(32) 



Remarkably, the cross section is identically zero when the direction of incidence is aligned 
with the field axis. Since, as will be discussed below, this semi-classical cross section 
accurately describes dipolar collisions in the temperature range currently accessible 
experimentally, there are a variety of observables which might test this angle- dependence 
of the total elastic cross section. While equilibrium properties of the gas would not be 
sensitive to the angle- dependence, non-equilibrium properties, such as transmission of 
fast dipoles through a trapped dipolar gas, would be expected to show strong dependence 
on the alignment of the beam with the field axis. 

Averaged over incident directions within a confined gas, one then expects 

(X e; 87T 87T 
D 2 = 



(33) 



3k 3KD 

which is the final result. As shown in Fig. 1, Glauber's method yields not only the correct 
scaling, but quantitatively reproduces the universal results discovered in close-coupling 
calculations pQ. 

The basic structure of this result was surmised by Gallagher [26J in a study of 
Rydberg-Rydberg collisions: from the uncertainty principle, an interaction of energy 
1/b 3 and lasting for a duration b/k, should satisfy l/b 2 k ~ 1 in scaled units, so that 
a ~ b 2 = 1/k. Using an isotropic — 1/r 3 interaction, DeMille [27] also applied the 
eikonal approximation and determined a cross section a/D 2 = 2n 2 /k. Kajita's Fourier 
technique [28] yields the high-velocity result a/D 2 = A0Tc^/2/3k. 
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3. The rise and fall of universal scattering 

Figure 1 has already presented the message of universality, in that the Born and eikonal 
limits are achieved in the appropriate energy ranges. However, the calculations in this 
figure were carefully selected to have zero scattering length, by choosing an appropriate 
cutoff radius r . By changing r , we are able to generate any scattering length a. 
Experimentally, the value of a can be altered by changing the electric or magnetic field 
strength. Changing ro could also change the scattering phase shift of any other partial 
wave. However, this is only likely for those partial waves whose centrifugal barrier lie 
below the collision energy, because these partial waves are the only ones to probe physics 
at the scale of r . 

The importance of the centrifugal barriers in determining the range of applicability 
of the low-energy (Born) scaling was emphasized in pEJ [5] , through adiabatic calculations 
which converge much more rapidly than close-coupling calculations at the lowest 
energies. The adiabatic curves have pronounced barriers (in all but the s-wave channel) 
separating repulsive centrifugal behavior at large-r from attractive dipolar behavior 
at small-r. The heights of the lowest barriers (approximately coincident with Ed) 
determine the range of energy over which the threshold behavior of the cross section is 
approximately constant. They also suggest the sensitivity to r at intermediate energies, 
where incident flux can surmount the barrier. 

Figure 2 shows the effect of changing r on the "universal" cross section from Fig. 
1. In Fig. 2a) is shown the result for even partial waves. The solid line is the a = 
result, and it amicably reaches the universal Born and eikonal limits. The other curves 
employ values of ro that produce scattering lengths of a = 0.1D (red) and a = —0.1D 
(blue). This change has made a significant difference in the low-energy limit, where now 
the Born approximation is merely a lower limit to the cross section. Notice also that 
for a < the cross section initially decreases with increasing energy, just as it does for 
alkali atoms [29] 

However, at higher energies E > Ed, this change in ro has no effect on a. One way 
to look at this is that the phase shifts have changed for many partial waves, but because 
there are so many of them added together to get the cross section, these changes average 
out. Another point of view is that the semiclassical scattering occurs at high impact 
parameter, and is thus indifferent to what happens at r = r . In the intermediate energy 
range, the change is still quite significant, since phase shifts are changing for only those 
few partial waves that skip over their centrifugal barriers. 

Figure 2b) shows the same circumstance, but for odd partial waves. The same three 
values of ro are employed here, and so the three curves are labeled by the scattering 
lengths from part a). In this case the cross sections always approach the Born limit, 
since there is no aberrant s-wave scattering to derail them. In the high-energy limit, 
too, a is again insensitive to r . It is only in the intermediate energy range that a 
small deviation is seen. This suggests that for odd partial waves the behavior of the 
cross section is indeed nearly universal. This would be true, for example, in collisions 
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Figure 2. Cross sections as in Fig. 1, but separated into contributions from even (a) 
or odd (b) partial waves. In each case, three different values of the cutoff radius tq 
are chosen, corresponding to three different s-wave scattering lengths a, as indicated in 
the legend. In each panel differences in the three curves demonstrate the breakdown 
of universality. Elastic scattering is thus much more universal for odd partial waves 
(identical fcrmions) than for even partial waves (identical bosons or distinguishable 
particles. 



Quasi- Universal Dipolar Scattering 



14 



of fermionic molecules (e.g., 40 K 87 Rb [3D]) in identical hyperfine states. 

To further emphasize the consequence of an s-wave scattering length, Figure 3 
reports the cross section a for the even partial waves as a function of Tq, in the threshold 
limit, using E/Ep = 10~ 3 . The range of r shown here corresponds to a complete cycle of 
the scattering length from zero, through infinity, and back to zero again. Consequently, 
the numerically evaluated cross section shows a resonance, at which point the cross 
section is determined by a 2 , not D 2 . Even away from the resonance peak, the s-wave 
contribution can significantly increase the cross section. We are thus led to conclude 
that universality at low-energies, in cases where s-wave scattering is allowed, is similar 
to the universality for atoms. Namely, the form of the cross section (independent of 
energy) is universal, but its value relies on a (field-dependent) scattering length that 
must be determined empirically. The Born approximation does provide a useful lower 
limit, however. For odd partial waves, while resonances exist, they are shape resonances, 
hence narrow at low energies and less likely to destroy universal behavior. 

The s-wave contribution to the Hamiltonian nominally vanishes, since Cqq — 0. 
However, it is not unreasonable that s-wave scattering has a strong influence near 
threshold. To see this, we evaluate an effective s-wave interaction at long range, via 
its coupling to the I = 2 partial wave, in second-order perturbation theory (compare 



where, in dipole units and using I = 2, the coefficient is C 4 = 4/3v5. This in turn leads 
to a characteristic length scale for the s-wave interaction, analogous to the characteristic 
van der Waals length, 



In dipole units, this is r 4 = 1.09-D, comparable to the dipole length itself. Based on this 
consideration, it is perhaps not too surprising that s-wave scattering plays a significant 
role. 

4. Dependence on incident angle 

Thus far we have focused on the cross section as integrated over all incident angles. One 
of the interesting aspects of the dipole-dipole interaction, however, is its anisotropy. 
The cross section cr to t(#i) may therefore depend on the angle 0$ = arccos(4 • S) of the 
incident collision axis, with respect to the polarization direction. This cross section is 
easily calculated numerically from (jSJ), and also from the useful eikonal estimate ( J32l) . 

To show the utility of the eikonal expression, we present in Figure 4 cr tot versus 
6i, for a collision energy E/Ep = 10 4 where the eikonal approximation should be 
fairly accurate. The total close-coupled cross section (black solid line) is the sum of 
contributions from even and odd partial waves. Note that the contributions from these 



Ref. j3U): 




(34) 




(35) 
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Figure 3. Cross section a for various values of the cutoff radius tq. The numerically 
determined cross section is always larger then the Born approximation, sometimes 
significantly higher, due to a large s-wave scattering length. 

two sets of partial waves are nearly equal here in the semiclassical limit where many 
partial waves contribute, and effects of dipole-indistinguishability are small. The angular 
distribution of each shows oscillations, but in the sum, representing distinguishable 
particles, the angular dependence is smooth. Moreover, for angles where the collision 
axis is orthogonal to the polarization axis, « 7r/2, the eikonal approximation (dotted 
line) is quite good. 

A major deviation occurs, however, for dipoles aligned parallel to the collision axis. 
Here the eikonal result calls for vanishing cross section, whereas the close coupling 
calculation yields a non-zero cross section. For = 0, the incident wave e lkl ' r is 
invariant under rotations about the field axis, and so contains only m = partial waves. 
It is not too surprising that a semi-classical analysis will break down for low-m states. 
Furthermore, for m = states, the wavefunction is large where the potential is strongest 
(near r = 0), so the eikonal assumption V/e << 1 is no longer valid. Most interesting 
about this deviation from eikonal behavior, is the importance of back-scattering from 
the strong potential in this geometry, as indicated by the pronounced difference between 
even and odd partial waves. Observations near Bi = will accordingly be most sensitive 
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Figure 4. Dependence of scattering cross sections <r tot on the incident angle 0i between 
the collision axis and the polarization axis of the dipoles. This calculation is performed 
at a high collision energy, E/Ejj — 10 4 . The distribution is well-approximated using 
the eikonal result except when the collision and polarization axes nearly coincide. 



to exchange scattering, and, presumably to short range physics. 

Cross section variations with the angle of incidence have also been studied at low 
energies [HE]. When the s-wave scattering length is negligible, a universal anisotropic 
distribution is obtained, entirely due to long-range scattering. However, when the s- 
wave scattering length dominates, near the peak in Figure 3, a completely isotropic 
distribution is found, as in the case of ultracold atomic collisions. Interestingly, this 
implies that effects of anisotropy are to be seen at the lowest temperatures only when 
the scattering length is small, and cross sections are accurately represented by the Born 
approximation. 



5. What this means for you 

Scaled units are fine for proving a theoretical point, as we have hoped to do here. 
However, since dipole length scales vary widely between different molecules and at 
different electric field strengths, it is also useful to consider specific examples that 
measure cross sections in cm 2 . 

Before presenting such an example, we first recapitulate our main results, cast in 
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Figure 5. Elastic cross section for scattering of pairs of fermionic °K 87 Rb molecules 
in identical internal states, averaged over incident directions. This calculation is based 
on the "universal" calculation that includes only dipole-dipole interactions. The top 
curve is the cross section for fully polarized molecules with dipole moment (fi) — 0.566 
Debye. The dipole is halved for each successively lower curve. 



terms of the explicit dimensionful factors. The Born result, valid in the ultracold limit, 



is 



MWW , , 2 

-Born = 1-11' ^4 ^ 47Ta 



-Born = 3.351 - . (36 ) 



K 

The semiclassical result, valid for cold collisions, E > Ed = fr 6 /M 3 (/ii) 2 (/i 2 ) 2 , is 




8tt (/ii)(/i 2 ) , .... 

<m = T ^—\l7^- (37) 

(Recall that our eikonal derivation does not distinguish between even and odd partial 
wave contributions. To a good approximation, both the even and odd contributions 
would be half this value.) These formulas show explicitly that the cross sections at 
low and high energy differ not only in their dependence on energy, but also in their 
dependence on the parameters - reduced mass and dipole moments - of the molecules. 

To give a concrete example, consider the ground-state 40 K 87 Rb molecules that were 
recently produced at temperatures of several hundred nK |30j. This molecule has a 
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dipole moment of 0.566 Debye, and, being a fermion, would collide only in odd partial 
waves if it is trapped in a single quantum state. We therefore plot in Fig. 5 the odd 
partial wave cross section computed above, but cast in realistic units for this molecule. 
The largest cross section corresponds to the full dipole moment, (//) = 0.566 Debye. 
Each successively lower curve divides the dipole moment in half from the previous one. 
For this reason, each low-energy cross section drops by a factor of 16 from the one 
above, while at high energy each cross section drops by a factor of 4. Because the 
dipole moment is something that can be changed by the application of a greater or 
lesser electric field, cross sections spanning this stunning range of magnitudes should be 
observable in experiments. 

Also interesting is the energy scale encompassed by this figure. In the experiment, 
the gas is trapped at a temperature of 350 nK. At low electric field values, hence low 
dipole moments, these molecules are in the ultracold regime, and scatter according to 
the Born prescription. At higher fields, however, Ejj approaches the temperature of the 
gas, and experiments might start to observe the non-universal behavior of the scattering. 

In summary, we have characterized the total scattering cross section for dipolar 
molecules, both in the cold limit Ed < E < B e , and in the ultracold limit E < Ed- The 
behavior of this scattering is universal for cold collisions, and nearly so for ultracold 
collisions. In the temperature regime intermediate between these two, universality 
breaks down, and empirical cross sections will likely reveal information about the 
intermolecular potential energy surface. 

The authors acknowledge the financial support from the NSF (J. L. B.), the ARC 
(C. T.), and the University of Kentucky (M. C.) 
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